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Flux flow properties for clean and dirty superconducting limits in strong pinning niobium 
films are studied. Measurements of electric field vs current density characteristics at high 
J values (J ~ l^fA/cm?) are successfully compared to theoretical models of flux flow and 
of flux flow instability in a large range of temperatures and magnetic fields. The non-linear 
regime at high dissipation is analysed in the frame of a modified Larkin Ovchinnikov model 
that takes into account a quasiparticles heating effect. This model elaborated by Bezuglyj 
and Shklovskij (BS) defines a transition magnetic field Bt above which the quasiparticles 
distribution became non uniform due to a finite heat removal from the substrate. From 
the BS model, we can deduce values of the non-equilibrium lifetime of quasiparticles r^p, 
which are 10 to 100 times shorter in the dirty sample compared to the clean one and whose 
temperature dependance is specific to the electronic nature of the Nb film. The study of the 
non-linear regime provides also a quantitative determination of the thermal transparency of 
the film-substrate interface. 

PACS numbers: 74.25.Fy ; 74.25. Qt ; 74.25. Sv ; 74.78.Db 

I. INTRODUCTION 

Magnetic vortices in the mixed state of a superconductor can move under the action of the 
Lorentz driving force Fl proportional to a transport current density J—Jc, where Jc is the pinning 
threshold. Vortex dynamics near Jc has been recently widely explored by direct methods such as 
magneto-optics^. Inversely, the regime at high vortex velocities is only deduced from electric 
field vs current density characteristics E{J) in accordance with the Josephsoi*^ equation E = 
—v^ X Ba, where Ba is the magnetic field inside the superconductor. A first experiment^ in 1964 
has interpreted a linear dependence between E and J as a steady fiux fiow (FF) of vortices. The 
FF regime is defined by the balance between Fl and a viscous drag force F^ = —rjppv^ where 
rjFp is the viscosity associated to vortex motion. A linear regime J = appE + Jc for which the 
viscous coefficient rjpF and the conductivity app {app = are constant is thus expected. 

The dissipation in a moderate "clean" superconductor at low magnetic fields and temperatures is 
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essentially due to currents crossing vortex cores^ : 

crpF ^ Bc2{T) 
aniT = 0)~ Ba 



k{T), (1) 



with (T„(T = 0) the conductivity of the normal state, Bc2{T) = and k(T) a temperature 
dependent positive constant representing the dissipation due to Cooper pairs around vortex cores. 
This equation can be simplified in: 

— = A{T)/Ba - k{T) (2) 

where A(T) is a temperature dependent coefficient. The expression of the conductivity apF for a 
"dirty" system adds the contribution of thermal effects^ in the vicinity of the vortex core : 

where is a positive function with a maxima at t = T/Tc ~ 0.5. 

At higher velocities v^, Larkin and Ovchinnikov'*' (LO) have calculated a non-linearity of the con- 
ductivity (Tff{vv) for temperatures T close to Tc and Ba <^ Bc2- High electric fields in the vicinity 
of the vortex core can change the electronic distribution and the electronic temperature Tgp leading 
to a decrease of r]pF when increases. E{J) characteristics become non- linear and end with a 
jump into the normal state before the depairing current density is reached. In the absence of ther- 
mal runaway, the flux flow instability occurs at a critical vortex velocity v* . To come to this critical 
point {J*,E*), the quasiparticles inside the driven vortex cores have gained enough energy from 
the electric field to escape from these normal regions and relax their energy into the condensate. 
This process is controlled by a non-equilibrium lifetime Tqp of the quasiparticles. This electronic 
leakage leads to a continuous shrinkage of the vortex core radius ^ and to a decrease of apF as a 
function of Vy. 

In the original LO theory, the non linear flux flow behavior is due only to the field induced 
change in the quasiparticle distribution function and not Joule heating: the system is in thermal 
equilibrium with the bath, characterized by a temperature Tq. Several experimental works in 
lo w^i^i^^ and higbii^i^ Tc superconductors have confirmed the validity of the LO model near the 
critical temperature, with however some discrepancies. 

It comes out that the sample heating during the dissipative flux flow is not negligible and can yields 
a thermal runaway^- before the occurrence of the FF instability. Bezuglyj and Shklovskij^^ (BS) 
have extended the LO theory in the thin film configuration by taking into account heating effects. 
In their model, the quasiparticles distribution function depends on the vortex density and on the 
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rate of heat removal from the film through the substrate. BS defines a transition magnetic field Bt, 
dependent on the thermal exchange coefficient h between the film and the substrate, under which 
the hypothesis of a uniform quasiparticles distribution (LO model) is still valid : 

= ^fi^, (4) 

where e is the electronic charge, ks the Boltzman constant and d the film thickness. For Ba » Bt, 
dissipation during the FF regime raises the electronic temperature Tqp and thermal effects govern 
the FF instability. This effect can be understood from a simple dynamical picture where the inter- 
vortex spacing becomes small enough to allow an influence of a vortex core on the other. In other 
words, the condensate keeps a memory in terms of quasiparticles energy (and temperature) from 
vortex passing. In contrast to the B-independent of the LO model, a v*{Ba) variation is now 
expectedii and takes the form : 

vl^hil-tf/^B:, (5) 

with n = —0.5. BS proposed a scaling law between critical parameters J* and E* for the heated 
quasiparticles at Tqp » Tq : 

- 4bt)){-r^r\ (6) 



^LO '-'LO 

where -E^q and J£q are critical parameters of the pure magnetic (LO) theory and z{bt) is a function 
of bt = BjBt. 

Recent work o^^i-^^ at low temperatures (T < 0.4Tc) suggest that thermal effects can diminish 
the superconducting order parameter and lead to an expansion of vortex cores rather than to a 
shrinkage. The quasiparticles heating reduces critical field Bc2{T) to Bc2{T*p) = Ba where the 
transition to the normal state occurs above v*. The Ba dependance of v* is the same than the one 
given by the BS theory. 

In this paper, we explore and compare the dynamic of a vortex lattice at high velocities in niobium 
micro-bridges for "clean" and "dirty" superconducting limits in a range of temperature above O.GTc. 
We reveal the influence of the electronic nature of superconductors on flux flow properties in the 
framework of the BS model. Values of h and the temperature behaviour of r^p, which are important 
intrinsic parameters for the development of electronic devices such as hot electrons bolometers, are 
also deduced. 



4 



TABLE I: Example of fundamental characteristics for "dirty" and "clean" Nb films 



Sample 


substrate 


Tc^^o (K) 


H,2(0) (mT) 




S.OBCS (nm) 


Ifree (nm) 


Nb dirty 




8.92 


4600 


12.12 


36 


3.19 


Nb clean 


AI2O3 


9.13 


1010 


0.59 


33 


65 


Nb dirty 


Si 


8.6 


4430 


9.9 


35 


3.9 



II. EXPERIMENTAL DETAILS 

Niobium thin films with thickness d ~ lOOnm are prepared by ion beam technique. Depositions 
are done either at ambient temperature (cool sample in dirty limit) or at 780 ° C (warm sample 
in clean limit) on Si and AI2O3 substrates. Films are protected by a silicon thin layer of 5nm 
thickness. Microbridges of 8.5 to lOfim width {w) are patterned to achieve a four points configura- 
tion measurement. Bridge lengths between voltage contacts are included between / = 800^m- and 
/ = 3mm. Pinning is strong in niobium films: values of critical current densities Jc are typically 
a few lO^A/cm^ at 0.7Tc. These sample parameters are within the range of expected values for 
Nb films. In all the experiments presented here, the magnetic field is applied perpendicular to the 
samples surface. More details are given in reference 17. 
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FIG. 1: E{J) curves at T = 7.SK (t ~ 0.9) for < Ba < 150toT (Bc2(T) = 420 ± 8 mT ) in a "dirty" 
niobium film. 



Current-voltage characteristics were measured through fast current sweeps {vj € [1; 250^4/5]). 
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An example is reported in Figure^ The electric field E and the current density J were determined 
from relations E = V/l and J = I/S where / and S = d x w are respectively length and section 
of microbridges. Parasitic thermal effects coming from contact resistance and classical Joule dis- 
sipation due to vortex motion are identified by performing experiments at different current sweep 
rates. To retain only non-equilibrium effects pertinent for the BS model, a high enough current 
rate vj where the parameters J* and E* are constant (see figure ISJ is applied. 

E(V/m) 
160 
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80 

40 



1.5 lO*" J*3 10*^ 4.5 10^ 6 10^ 
J (A/cm^) 

FIG. 2: Variations of J* and E* vs vi at T = 6K and Ba = 20mT in a "dirty" niobium film. See Fig. 1 
for the experimental determination of J* and E*. These critical parameters are independent of the sweep 
rate above 20A/s. 




III. RESULTS AND DISCUSSION 

Figure ^ depicts a serie of typical {E, J) curves in a "dirty" niobium sample at t ~ 0.9 for 
< b{t) = Ba/Bc2(T) < 0.35. For current densities J > Jc, a linear regime E (x J occurs, 
which corresponds to a constant viscosity 7]pp (see also Fig. 3a). At higher currents, the response 
becomes non- linear (decreasing rjFp) and finally ends by a flux flow jump at low magnetic fields. 
This general behavior is reported for both kinds of superconducting films. 

The flux flow can be characterised by two thresholds of vortex speed, one defining the onset (vmin) 
of the linear FF regime, the other its end (vmax)- It is interesting to note that these FF velocity 
thresholds are independent of the applied field, at least in the intermediate field range 10mT-70mT 
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as shown in Fig. 3b. A similar behaviour is found in the clean limit film with however much lower 
Vmax (about 110 m/s). 
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FIG. 3: Right : E(^J) curves at T = T.Sif {t ~ 0.9) for Ba = 30, 40, TOtoT in a "dirty" niobium film. Left : 
Magnetic field dependence of the thresholds of the vortex velocity Vmin and Vmax ■ 



The constant value of the onset of the linear regime can be interpreted by a dynamic phase 
transitio n^^i^^ under current from the motion of the amorphous vortex configuration at J < Jc(u = 
Vmin) (plastic flow) to the motion of the vortex crystal (moving crystal) at J > Jdv = Vmin)- The 
crystalliszation current is expected to exceed depinning current J^ep and is related to lattice defect 
concentration. Alternatively, a transition from vortex creep to fiux flow with increasing Lorentz 
force can be considered. The existence and the value of the maximum vortex speed ending the 
linear FF regime will be discussed below in the section dedicated to the non-linear behaviour at 
high dissipation. 

The conductivity app in the linear regime is treated within the ranges 0.02 < b < 0.5 and 
0.65 < t < 0.95 in the framework of the fiux flow model. Figure 0] analyses the Ba dependence of 
the ratio cpp/an and shows the good accordance of experimental data with the A{T)/Ba — k{T) 
formula for both kind of samples. For "clean" systems, data on almost two decades of dissipation 
agree with equation (2), giving A(T) absolute values shifted by only 3.5% from the Hc2{T) curve 
deduced from AC resistivity measurements (see inset of figure 0J). This result validates the approx- 
imation of a normal core with a radius of the order of £,{T) inside which the dissipation occurs. 
On the contrary, the coefficient A{T) for "dirty" films doesn't not follow the temperature variation 
given by equation (2) and displays a maximum at t ~ 0.85. For temperatures close to Tc, A(T) 
decreases toward the normal conductivity but remains above the values, meaning that the ratio 
app/uN is higher than expected, i.e. the dissipation inside the core is lower than what is found 
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in clean samples. Larkin and Ovcliinniko\»^ have calculated app in the theoretical frame of the 
Ginzburg-Landau theory. Their formula depending on T and Ba doesn't fit our present data for 
"dirty" samples. A complementary theoretical work in thin film configuration where the thickness 
is of the same order than the magnetic penetration length A is needed. At high magnetic fields, a 
recent theoretical work^ has successfully fitted the experimental curves of Figure 1 with a modified 
time-dependent Ginzburg Landau equation. In summary, dissipative phenomena governing the FF 
regime are very different between "clean" and "dirty" films although a clear linear E{J) regime is 
observed in both cases. 

Looking now at higher dissipation, properties of the jump into the normal state are discussed 




FIG. 4: Ratio (JpF/crn as function of b{t) at i ~ 0.9. Sample (1) and (2) are respectively "clean" and 
"dirty" films on saphire substrate, and sample (3) is a "dirty" film on a Si substrate. Dashed lines represent 
fits to equation A{t)/Ba ~ k{t) with A{t) and k{t) as fitting parameters. Inset : Comparison of A{t) points 
and Hc2{t) curves (solid lines). For sample (1), Hc2{t) is corrected by a factor 1.035. 



within the BS model for which the finite thermal exchange h between the superconducting film and 
the substrate plays a role. Values of critical coordinates {J*,E*) are confronted with the scaling 
law defined by equation (6). First, the field Bt is directly deduced from the z{bt) dependence of the 
dissipative power P* = E* x J* . Inset of Figure [HI shows good accordance between experimental 
data and theoretical curves with two ajustable parameters Bt and P£q = E^q x J£q. In a second 
step, E^Q and J^q are extracted from the Bt (i.e. z{bt)) dependance of J* and E* according 
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to the BS model. The curves J*{z) and E*{z) fohow the dependence predicted by equations 33 
and 35 of referenceii with only one adjustable parameter per equation, respectively J^q and Ef^Q. 
This analysis is performed in a large range of T and Ba as seen in Figure El which displays the 
expected scaling law predicted by equation ^ The extracted values of Bt are about 5 — lOmT 
in our experimental range for all samples, which corresponds to vortex densities around 2 to 5 
vortex / fim^ . These low densities tend to show that the regime of quasiparticles heating occurs 
in a large field domain even for low temperatures. The field dependence of v* in the form B^ in 




J7J*L0 

FIG. 5: Comparison between {E* , J*) data and universal scaling law y = 1/x. Open symbols are for the 
"dirty" film at T = 6; 6.5; 7; 7.5; 8K. Triangle symbols are for the "clean" sample at T 6.5; 7K. Inset : 
P* — -P£o[l — ^(^)] variation as a function of Ba a,t T — 6.5K for the "dirty" film. The two parameters 
P£q and Bt are deduced from this fit. 

both superconducting limits (Figure EJ, which is not predicted in the frame of the LO model, gives 
another proof of the quasiparticles heating mechanism in our samples. The experimental n values 
are close to 0.4, which is in a reasonable agreement with the BS model, although some refinements 
of the theory are needed here. The highest velocities v* are found in the dirty limit. For "dirty" 
films, the theoretical relation fits exactly the experimental v*(T) variation (Inset of Figure IHl) in the 
range B^ > Bt by taking a constant heat exchange coefficient h versus temperature. In contrast, 
v*{T) for "clean" niobium shows that h will vary a lot with temperature. 

Taking into account the good accordance between experimental data and theory, Tgp{T) is ex- 
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tracted from the product P^q = J^q x E*j^q and the normal conductivity o"„ (equations 34 and 
36 of referencei^ giving r^p = 2.67-^^^/c_bTc(1 — t)). This Tqp calculation has the advantage to 
not depend on Fermi velocity vp which is difficult to know precisely. Two very different micro- 
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Ba (mT) 

FIG. 6: Extracted v*{Ba) values for "clean" and "dirty' films at T = 7.5K. Solid fines fit data with B" : 
n ~ 0.35 and n ~ 0.38 respectively for "clean" and "dirty" film. Inset : corresponding v*{T) for Ba ~ 20mT. 
Solid line represents (1 — T/Tc)^-^^ fit curve. 

scopic mechanics are revealed with the rgp(T)dependence for both superconducting limits. For 
"dirty" niobium films, quasiparticles lifetime doesn't change much with temperature and lays in 
the range 30 — 50 ps, this order of magnitude being in accordance with data obtained from magneto- 
conductance measurements^^ and photonic Nb detectors^^ (figur^Tj). 

This order of magnitude is also compatible with a simple calculation considering that vortex mo- 
tion is characterised by pair breaking upstream from the core and by quasiparticles recombinaison 
downstream. Quasiparticles inside the core can relax their extra energy given by the electric field 
only if their characteristic time of presence in this normal region is larger than Tqp. When the 
vortex motion becomes too fast, this energy relaxation can not occur anymore and an effective 
increase in quasiparticle energy that will ultimately lead to the instability, takes place. In this 
description, Tqp obeys the simple law S,/vmax = Tqp. As ^ ?a lOnm in the dirty sample at 7.8K, 
we obtain with Vmax = 480m/s (see Figure 3b), a characteristic time Tqp = 6.10^^^s, in perfect 
agreement with the value obtained from the BS model. 
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In the clean sample the same analysis yields Tqp = 5.10" s for this reduced temperature t=0.9, 
which is again perfectly compatible with our experimental data {(,t=7.5K = 58nm, Vmax = llOm/s). 
The quasiparticle energy relaxation process is however ruled by different mechanisms in the two 
samples. When the mean free path is shorter than the coherence length (dirty sample), the energy 
relation occurs within the vortex core. In the opposite situation (clean sample), the quasiparticles 
can relax their energy given by the electric field only after covering a distance within the core 
equivalent to several i^. In the microscopic picture where the quasiparticle energy rise comes from 
Andreev reflections at the vortex core walls, one can see that an increase of the quasiparticles 
energy can occur below Vmax in clean samples but it is only above Vmax that quasiparticles can 
escape from the core, leading to its shrinkage. This efficient process of quasiparticles energy rise, 
even at low in clean samples, is coherent with the fact that Vmax is here lower than in dirty 
systems. 

The above discussion has to be linked to the study of the Tqp{T) dependance which is related to 
the physical mechanisms attached to this parameter. For example, quasiparticles recombination to 
Cooper pairs with emission of phonons follows a law in e^^'<'^)/^BT with the temperature variation 
of the superconductive gap being A(T) ~ A(0)(1 — t)^-^. Figure 7a shows that the Tqp{T) variation 
in the clean sample conveniently fits with a recombination model^. The only adjustable parameter 
A(0) is found to be about IksTc, higher than the BCS value A(0) ~ l.TGfc^Tc. Such a deviation 
from the weak coupling BCS limit has already been reported in the case of Nb thin films^^. On 
the contrary, a quasi-constant Tgp(T) variation in dirty samples show that scattering mechanisms 
dominate and that quasiparticles are ejected from the core at energies 5E » A(0), as suggested 
by S.B. Kaplan et ali^. The predominance of scattering mechanisms is coherent with the fact that 
quasiparticles relax their energy in the core on a distance I free < ^- The limiting mechanism of 
the instability is thus the inelastic scattering within the core for the dirty sample. Recombination 
process on the other hand describes well the fact that quasiparticles scattering within the core is 
not the predominant phenomena which restrains the electronic leakage. The limiting mechanism 
will take place outside the core, which is coherent with a recombination mechanism. Although 
an extensive microscopic calculation would be needed to describe the physics of a vortex core 
in motion, we can already see that the phenomenology of clean and dirty samples will be quite 
different and is indeed associated to distinct behaviours. 

Knowing Tqp and Bt, heat transfer coefficients h are extracted from equation (4). Values of h are 
of the same order of magnitude (around 2.7 to 4.7 W/cm'^K) for thin films on AI2O3 substrates for 
both superconducting limits. These values are one to two order of magnitude lower than h values for 
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FIG. 7: Variations of duration Tqp{T) (Fig. a) and coefficient h{T) (Fig. b). Dotted line in Figure a 
fits experimental data with recombination model e^^^'^^^'^^^. Squares symbols are data for Au/Nb/ Au on 
saphire substrate. 

epitaxial BiSrCaCuO films on SrTiOs substrates^^ deduced from the first experimental analysis 
based on BS model, or for epitaxial YBaCuO films on saphire substrates^I, but are 10 to 100 times 
higher h found for YBaCuO monocristals^'^ only glued. Present values seem realistic although a 
higher heat transfer coefficient would have been expected in the case of the clean (oriented) thin 
film. We however find that h is independent on the electronic nature of Nb films and shows high 
values of thermal resistance between substrate and superconducting film. In this way, the thermal 
contact is improved with the insertion of a metallic layer between substrate and Nb film (Fig. 7b). 
As a final remark, we can compare our results obtained in the present work to those reported in 
our previous paper^ considering only a non thermal LO process. A significant differences is found 
for the clean samples where the Tqp{T) variation extracted from the BS model now follows a law 
relative to a recombination process with a very reasonable A(0). 
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IV. CONCLUSION 

We have compared flux flow properties in clean and dirty superconducting limits of niobium thin 
films. For "clean" films, dissipation is essentially due to currents crossing vortex core with a radius 
^(T). On the contrary, FF properties is governed by thermal effects in the vicinity of the vortex 
core in "dirty" films. We have identified a non-linear regime at high dissipation in the frame of a 
modified Larkin Ovchinnikov model that takes into account a quasiparticles heating effect. Non- 
equilibrium lifetime of quasiparticles Tqp in "clean" films is interpreted as a recombination process. 
Values of Tqp are about 10 to 100 times (few lOps) shorter in "dirty" niobium films and reveal 
different microscopic mechanism for quasiparticles relaxation suggesting high performances for hot 
electrons bolometers. We hope that this study will motivate future experimental and theoretical 
works to understand in details FF properties and non-linear regime in "dirty" superconducting 
limit of thin films. 
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